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Abstract 



In this paper we express the velocity dependent, spin dependent heavy quark 
potential V q q in QCD in terms of a Wilson Loop W(r) determined by pure Yang 
Mills theory. We use an effective dual theory of long-distance Yang Mills theory 
to calculate W(T) for large loops; i.e. for loops of size R > Rft- (Rft is 
the flux tube radius, fixed by the value of the Higgs (monopole) mass of the 
dual theory, which is a concrete realization of the Mandelstam 't Hooft dual 
superconductor mechanism of confinement). We replace W(r) by W e ff(T), given 
by a functional integral over the dual variables, which for R > Rft can be 
evaluated by a semiclassical expansion, since the dual theory is weakly coupled 
at these distances. The classical approximation gives the leading contribution to 
W e ff(T) and yields a velocity dependent heavy quark potential which for large R 
becomes linear in R, and which for small R approaches lowest order perturbative 
QCD. This latter fact means that these results should remain applicable down to 
distances where radiative corrections giving rise to a running coupling constant 
become important. The spin dependence of the potential reflects the vector 
coupling of the quarks at long range as well as at short range. The methods 
developed here should be applicable to any realization of the dual superconductor 
mechanism. They give an expression determining W e ff(T) independent of the 
classical approximation, but semi classical corrections due to fluctuations of the 
flux tube are not worked out in this paper. Taking these into account should 
lead to an effective string theory free from the conformal anomaly. 
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1 Introduction 



In this paper we give expressions for the heavy quark potential in QCD using an effective dual 
theory of long distance Yang-Mills theory. This work goes beyond a previous treatment M 
where the quark motion was treated semi-classically and where the dual theory was consid- 
ered only at the classical level, and provides an independent approach to the problem of the 
heavy quark potential. 

In Section two we give the formulae for the heavy quark spin dependent velocity depen- 
dent potential V q q obtained in refs. t 2 ' 3,4 ' in terms of a Wilson loop W(T). This expression 
extends previous work of Eichten and Feinberg,' 5 ' Peskint 6 !, Gromes and others'^ to include 
the velocity dependent spin independent part of the potential. The problem of the heavy 
quark potential is then reduced to the problem of calculating W(T) in pure Yang-Mills 
theory. All momenta, spins, masses and quantum mechanical properties of the quarks ap- 
pear explicitly in the formulae! 2 ! relating V qq - to W(Y). The size of the loop V fixed by 
the classical trajectories of the moving quark- ant iquark pair provides a length scale R (the 
quark- ant iquark separation) and we use the dual theory to evaluate W(T) for R > Rft, the 
radius of the flux tube that forms between the moving quark antiquark pair. 

In Section three we describe the dual theory and then show how to calculate the Wilson 
loops of Yang Mills theory at long distances (large loops). This is done by replacing W(T) by 
W e ff(T), a functional integral over dual potentials which are the fundamental variables 
of the dual theory. We then obtain the spin independent part of the heavy quark potential 
directly in terms of W e ff(T). Finally we discuss the relation of the dual theory to recent 
work' 7,8 ' on the use of electric-magnetic duality to determine the long distance behavior of 
certain supersymmetric non-Abelian gauge theories. 

In Section four we give explicit expressions for the spin dependent part of the heavy 
quark potential in terms of quantities determined by the dual theory. Since the theory is 
weakly coupled at large distances, W e ff(T) and hence V qq can be evaluated by a semiclassical 
expansion. 
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In Section five we calculate W e ff(T) in the classical approximation to the dual theory. We 
show how these results yield the dual superconducting picture of confinement and discuss 
their relation to the "modified area law" model for W(T), proposed in ref.' 9 '. Finally we 
remark how recent progress on quantization around classical vortex solutions' 10 ' may be 
useful for calculating corrections to W e ff(T) accounting for fluctuations of the length of the 
flux tube. 

In Section six we use the results of Section four and the classical solution to the dual 
theory to obtain the spin dependent part of the heavy quark potential. This calculation 
gives a contribution to V qq not obtained previously' 11 ] and yields a simplified expression for 
the spin orbit potential which reflects the vector nature of both the short range force and 
the confinement force. 

In the conclusion we point out that the results presented here should be regarded more 
as consequences of the dual superconductor picture in general rather than of our particular 
realization of it.' 12 ' 

2 The Heavy Quark Potential in QCD 

To obtain the heavy quark potential V qq we' 2 ' make a Foldy Wouthuysen transformation on 
the quark-antiquark Green's function and show that the result can be written as a Feyn- 
man path integral over particle and anti-particle coordinates and momenta of a Lagrangian 
depending only upon the spin, coordinates, and momenta of the quark and antiquark. Sepa- 
rating off the kinetic terms from this Lagrangian one can identify what remains as the heavy 
quark potential V qq . (Closed loops of light quark pairs and annihilation contributions were 
not included.) 

The terms in V qq of order (quark mass) -2 are of two types; velocity dependent Vvd an d 
spin dependent Vsd- The full potential V qq is then 

V qq = Vo(R) + V VD + V SD , (2.1) 
where V (R) is the static potential. These potentials are all expressed in terms of a Wilson 
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Loop VT(r) determined by pure Yang-Mills theory, given by 

_ J VAe iS ™( A krPexy(-ie & dx^A^x)) 

vv y L )- jvAe iS YM(A) ■ 

The closed loop T is defined by quark (anti-quark) trajectories zi(t)(z 2 {t)) running from 
yi to Xi(x2 to y 2 ) as t varies from the initial time U to the final time tf. The quark (anti- 
quark) trajectories z\(t)(z 2 (t)) define world lines ri(r 2 ) running from to tf(tf to t{). The 
world lines T\ and along with two straight lines at fixed time connecting y\ to y 2 and 
x\ to x 2 , then make up the contour F (see Fig.l). As usual A^(x) = ^XaA^x), tr means 
the trace over color indices, P prescribes the ordering of the color matrices according to the 
direction fixed on the loop and Sym(A) is the Yang-Mills action including a gauge fixing 
term. We have denoted the Yang-Mills coupling constant by e, i.e., 

a s = (2.3) 

The spin independent part of the potential, Vq + Vvd, is obtained from the zero order 
and the quadratic terms in the expansion of ilogW(T) for small velocities z\(t) and z 2 (t). 
This expansion has the form: 



ziogty(r) = dt U(R(t)) + E E ^mf{Rm-it)) ■ (2.4) 

where R(t) = Zi(t) — z 2 (t), and 



i,j=l k,£=l 



v VD =i: e ztm k /(Rm-(t). (2.5) 

i,j=l k,l=l 

(i\ogW{T) has an expansion of the form ()2.4j) only to second order in the velocities.) The 
expression (|2.5jl for Vvd follows from the same argument used to identify Vq(R) as the 
velocity independent term in the expansion ()2.4jl . We can write eq. ()2.4j) in the form 

i\ogW(T) = - f 1 dtL^z^z,,^), (2.6) 
Jti 

where 

— Li — V (R) + V VD (2.7) 



is an effective interaction Lagrangian for classical particles moving along trajectories <?i(t) 
and z\(t) with gauge couplings e(— e) and we can then interpret ilogVTfT) as an effective 
action describing the motion of classical particles after elimination of the Yang-Mills field. 

The spin dependent potential Vsd contains structures for each quark analogous to those 
obtained by making a Foldy Wouthuysen transformation on the Dirac equation in an external 
field F^ XT , along with an additional term Vss having the structure of a spin-spin interaction. 
We can then write 

VsD — ^LS + ^Thomas + Vparwin + Vgs, (2.8) 

using a notation which indicates the physical significance of the individual terms (MAG de- 
notes magnetic). The first two terms in eq. ()2.8|) can be obtained by making the replacement 

F^ XT (x) - «i^(*)», (2.9) 

in the corresponding expression for the interaction of a Dirac particle in an external field, 
where 

/{f(A))) - SVAe i ^"WtrP{e X j>[-ief r dx>>A ll (x)]f(A)} 

UJl >'' SVAe^M^HrPe^-iehdx^A^x)} ' 1 ' ] 

and 

F^ u = d»A v - d u A^ - ie[A^ A v ], (2.11) 

i.e. ((F^v^x))) is the expectation of the Yang-Mills field tensor in the presence of a quark 
and anti-quark moving along classical trajectories z\(t) and z\(ti) respectively. 
The explicit expressions for V£g AG and Vrhomas obtained in ref.^ arC 

/ dtVt s AG = £ — / dWSjaPUx))), (2-12) 



j=l 3 



3 Here and in the following J r dx tI f tl (x) = (— J t * dt(f (zj) — ij ■ f{zj)), where Zj — (t,Zj(t)). 
The factor (— accounts for the fact that world line runs from tf to fy. We also use the notation 

2£ = (f,3(f)). 
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and 
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/ dtV Thomas = -E A / dx^Sy^dF^x))), (2.13) 
where 

Sj is the spin matrix, and rrij is the mass of the jth quark. Because the expression for 
Vrhomas contains an explicit factor of l/m?, the integral over the trajectory of the jth quark 
f T dx^{{F fJlr (x))} can be replaced by (—iy +1 f t fdt ({F 0r (Zj)}} evaluated for static quarks. 
This gives the usual expression for Vrhomas in terms of the derivative of the central potential 
(see Section four). The expression for V^| AG on the other hand contains only a single power 
of ^- and J r . dx a ((p£ a (x))) must be evaluated to first order in the quark velocities. There 
results the usual magnetic interaction of the spin of the jth quark with the expectation value 

«^;)»- 

The expression for Vr) arwin is 

2 

/ dtV Darwin = -E A / dx»((D v F^(x)}), (2.15) 

where 

D v F Vik = d v F Vji - ie[A\ F ni ]. (2.16) 
Again because of the explicit factor of -\, the integral over the trajectory Tj of the jth 

i 

quark is evaluated for static quarks. 

The final term Vss i n eq. ()2.8|) is given by 

/ Vssdt = -i E -^— T * I dx» I dx'°S*S* 
J 2 .fy , hi jin y J\\ Jr., J J 

(((F E „(x)F ka (x'))) - m,(x)))((F ka (x')))), (2.17) 

where T s is the spin time ordering operator along the paths I\ and I^, and the averages 
are evaluated for static quarks. The terms j ^ f in eq. (|2.17|) give a spin-spin interaction 



proportional to l/m 1 m 2 while the terms j = f in eq. ()2.17|) give a spin independent term 
proportional to + The spin ordering is relevant only for these latter terms. 

We have thus obtained the explicit expression (2.8) for the spin dependent potential as 
a sum of terms depending upon the quark and antiquark spins, masses and momenta with 
coefficients which are expectation values (( )) of operators computed in Yang-Mills theory in 
presence of classical sources generated by the moving quark-antiquark pair. We now show 
that these expectation values can be obtained as functional derivatives of i\ogW(T) with 
respect to the path, i.e., with respect to the trajectories Zi(t) or Z2(t). For example consider 
the change in W (T) induced by letting 

zi{t) —* Zi{t) + Sz\(t), where 5z\{ti) = 5z\(tf) = 0. 

Then from the definitions (2.2) and (2.10), it follows that 

5zlog W(T) = - e jT' 5 -^^((F»„(zi)}), (2-18) 

where 

SS^izx) = (dzfSzt - dz»5z?). (2.19) 

Eq. (g3BD then gives 
and similarly one can get 

6ilogW(T) 
e((F^)}} = SSlv{zi) ■ 

Varying the path z 2 {t) in eq. ()2.2()j) gives 

e 2 ({(F, u ( Zl )F p(7 (z 2 ))) - ((F, u ( Zl ))}((F pa (z 2 )}}) = ie—J^—^F^fa))). 

db ^ (2.21) 

The first and second variational derivatives of VK(T) then determine the expectation val- 
ues of F pu needed to evaluate Vsd- Furthermore, we show in an appendix that ((D u F u ^(x))) 
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appearing in V£, arw i n can also be expressed in terms of variational derivatives of VK(r). The 
Wilson loop VT(r) which is determined by pure Yang-Mills theory then fixes the complete 
heavy quark potential V oS . Thus, up to order f r ) the dynamics of a quark anti- 

J ^ r yy \ quark mass J J ^ 

quark pair in QCD is completely fixed by the dynamics of Yang-Mills theory. The properties 
of the quark spins, masses, etc., appear only as given kinematic factors in the terms defining 
the heavy quark potential. 

The result (2.8) for Vsd is a consequence of the vector nature of the QCD interaction 
and contains precisely the same dependence upon the quarks spins, masses, and momenta as 
in QED. For example in eqs. (|2.12j) and (j2.13j) there is the usual vector coupling of quarks 
to ((F flv (x)J). The long (short) range part of Vsd is determined by the behavior of this field 
at long (short) distances. Both have the same vector coupling. 

This expansion as it stands is applicable only to calculating the potential between heavy 
quarks. The essence of the constituent quark model is that the same potential can also be 
used to calculate the energy levels of mesons containing light quarks with constituent masses 
fixed by hadron spectroscopy. The assumption is that the principal effect of the light quark 
dynamics can be accounted for by giving the light quarks effective masses which become the 
parameters of the constituent quark model. 

Finally we note the following " modified area law" proposed in ref. ^ : i log W(r) is written 
as the sum of a short range (SR) contribution and a long range (LR) one: 

i\ogW(T) =i\ogW SR (T) +i\ogW LR (T), (2.22) 

with i log W SR (T) given by ordinary perturbation theory and 

t\ogW LR (T) = aS mm} (2.23) 

where S m i n is the minimal surface enclosed by the loop T and a is the string tension. We 
will see in Section five the relation of this ansatz to the predictions of the dual theory. 
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3 The Dual Description of Long Distance Yang-Mills 
Theory 

The dual theory is an effective theory of long distance Yang-Mills theory described by a La- 
grangian density C e ff in which the fundamental variables are an octet of dual potentials C M 
coupled minimally to three octets of scalar Higgs fields Bi carrying magnetic color charge. 
(The gauge coupling constant of dual theory g — ^ where e is the Yang-Mills coupling 
constant.) The monopole fields Bi develop non- vanishing vacuum expectation values B 0i 
(monopole condensation) which give rise to massive C M and consequently to a dual Meiss- 
ner effect. Dual potentials couple to electric color charge like ordinary potentials couple to 
monopoles. The potentials C M thus couple to a quark anti-quark pair via a Dirac string con- 
necting the pair. The dual Meissner effect prevents the electric color flux from spreading out 
as the distance R between the quark anti-quark pair increases. As a result a linear potential 
develops which confines the quarks in hadrons. The dual theory then provides a concrete 
realization of the Mandelstam 't Hooft' 13 ! dual superconductor picture of confinement. 

Because the quanta of the potentials are massive, the dual theory is weakly coupled 
at distances R > (M being either the mass of the dual gluon or of the monopole field) 
and a semi-classical expansion can be used to calculate the heavy quark potential at those 
distances. The classical approximation gives the leading contribution to functional integrals 
defined by £ e ff, m contrast to the functional integrals of Yang-Mills theory where no single 
configuration of gauge potentials dominates W(T). The duality assumption that the long 
distance physics of Yang-Mills theory depending upon strongly coupled gauge potentials 
is the same as the long distance physics of the dual theory describing the interactions of 
weakly coupled dual potentials C M and monopole fields Bi forms the basis of the work of this 
paper. 

Before writing down the explicit form of £ e ff, we first show how to calculate W(T) for 
Abelian Gauge theory using the dual description of electrodynamics^ 14 ] , which describes the 
same physics as the original description at all distances. We consider a pair of particles with 
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charges e(— e) moving along trajectories i' 1 (t)(i* 2 (t)) in a relativistic medium having dielectric 
constant e. The current density j^(x) then has the form 

f{x) = e£dz tl 6(x-z), (3.1) 

where T is the world line described in fig.l. In the usual (electric) description this system 
is described by a Lagrangian 

CaU) = -\{d a A p - d p A a f - j a A a , (3.2) 

from which one obtains the usual Maxwell equations 

d a t{d a A p -dpA a ) =j p . (3.3) 

If the (wave number dependent) dielectric constant e —>■ at long distances, then we see 
from eq. (3.3) that A^ is strongly coupled at long distances (anti-screening). From (j3.1|) and 
()3.2j) we have 



dxC A (j) = - / dx A v ^ - e i dzPA^z). (3.4) 



t{dpA v - d v A v - 

4 " Jt 

The functional integral defining iy(r) in Abelian gauge theory 



I'D A e i I dx l c AU)+£-GF] 
W{T) = j VA ^Jd x[ c A u=o) + c GF V (3 ' 5) 

where Cqf is a gauge fixing term, is gaussian and has the value 

W(T) = e^Ir^Ir / ^ 2 ^, (3.6) 

where D^ v is the free photon propagator and where self energies have been subtracted. 
Because of current conservation the result ()3.6|) is independent of the choice of gauge. Let- 
ting e = 1 and expanding i\ogW{T) to second order in the velocities, as in eq. ()2.4|) . 
gives the Darwin Lagrangian Lp describing the interaction of a pair of oppositely charged 
particles 1151 ' 1161 

e 2 1 e 2 



J D 



AttR 2 AnR 



.. , (vi-R)(v 2 -R) 
Vl -v 2 + 



(3.7) 



In the dual description we consider first the inhomogeneous Maxwell equations, which 
we write in the form: 



-dP5&f^=j m (3.8) 

where is the dual field tensor composed of the electric displacement vector D and the 
magnetic field vector H: 

Gok = Hki Gem = £imnD n . (3.9) 

Next we express the charged particle current in eq. ()3.8|) as the divergence of a polarization 
tensor G^ u , the Dirac string tensor, representing a moving line of polarization running from 
the negatively charged to the positively charged particle, namely' 14 ' 

/f dy a Dy 
da J dT——5(x-y(a,T)), (3.10) 

where y a (a, r) is a world sheet with boundary T swept out by the Dirac string. Then M 



and the solution of eq. (3.8) is 



G^v — d^Cy — d u C^ + G^, (3.12) 



which defines the magnetic variables (the dual potentials C M ). (With eqs. (3.11) and (|3.12j) 
the inhomogeneous Maxwell equations become Bianchi identities.) 

The homogeneous Maxwell equations for E and B, which we write in the form 

d a (fiG afS ) = 0, (3.13) 

where \i — - is the magnetic susceptibility, now become dynamical equations for the dual 
potentials. These equations can be obtained by varying in the Lagrangian 

£c{G%) = -~/iG Mi ,G"-, (3.14) 
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where G^ u is given by eq. (|3.12|) . This Lagrangian provides the dual (magnetic) description 
of the Maxwell theory (3.2). In the dual description the Wilson loop W(r) is given by 

W(T) = J " , , a -. (3.15) 

Evaluating the functional integral ()3.15j) by completing the square gives 

W{T) = e-*!^W**<>*\ (3.16) 

where G al3 (x) is the dual field tensor ()3.12|) with = determined from the solution of 
eq. (|3.13|) . which has the explicit form 

d a fx(d a C^ - dpC?) = -d a fxG%. (3.17) 

Inserting 

G aP {x) = -\e a ^{d x f dyV^{x - y)j^y) - d a [ dyV w (x - y)jp(y)) 

1 J J (3.18) 

into 1)3. 16|) . integrating by parts, and using eq. (3.11), we obtain the same result (eq. ()3.6|) 
with - — > fx) for the Wilson loop (3.15) defined in the magnetic description as we had 
obtained for the Wilson loop defined in the electric description. We then have two equiva- 
lent descriptions at all distances of the electromagnetic interaction of two charged particles. 
(Note, however, that if e — > at long distances, then \i — > - — > oo and the dual potentials 
determined from eq. ()3.13|) are not strongly coupled at long distances unlike the potentials 
v4 M determined from eq. (|3.3jh ) 

We now return to C e ff, which in absence of quark sources has the form' 1 ' 

1 1 



■ i G^G„ + -(^B J ) 2 



W(Bi), (3.19) 



Ceff = 2tr 
where 

V ti B i = d^B i -ig[C li} B i ], (3.20) 

G, u = d,C u - d u C, - ig[G» C u ], (3.21) 
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2tt 

9=—, 
e 

Cft and Bj are SU(3) matrices, and W(Bi) is the Higgs potential which has a minimum at 
non-zero values B j which have the color structure 

Boi = B \ 7 , B 2 = B (— A5), B 3 = -80^2- (3.22) 

The three matrices A7, — A5 and A2 transform as a j = 1 irreducible representation of an 
SU(2) subgroup of SU(3) and as there is no SU(3) transformation which leaves all three 
B i invariant the dual SU (3) gauge symmetry is completely broken and the eight Goldstone 
bosons become the longitudinal components of the now massive C M . 

The basic manifestation of the dual superconducting properties of £ e // is that it gen- 
erates classical equations of motion having solutions' 1 ^ carrying a unit of Z% flux confined 
in a narrow tube along the z axis (corresponding to having quark sources at z — ±00). 
(These solutions are dual to Abrikosov-Nielsen-Olesen magnetic vortex solutions' 18 ' in a su- 
perconductor). We briefly describe these classical solutions here in order to specify the color 
structures that enter into the subsequent treatment of the dual theory with quark sources 
which is not restricted to the classical approximation. We choose a gauge where the dual 
potential is proportional to the hypercharge matrix Y — 

C M = C^Y. (3.23) 

As a consequence the non-Abelian terms in the expression (j3.21|) for the dual field tensor 
G^j, vanish. 

We choose Higgs Fields Bj having the following color structure: 

Bx = J B 1 (a;)A 7 + J B 1 (x)(-A 6 ) 
B 2 = B 2 {x){-\ 5 ) + B 2 {x)\ 4 

B 3 = B 3 (x)\ 2 + B 3 (x)(-\ 1 ). (3.24) 
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With this ansatz the Higgs potential W turns out to be 

+ 7(|0 1 | 2 + |0 2 | 2 + |0 3 | 2 -3 J B 2 ) 2 } , 

(3.25) 

where 

<j> i {x) = B i {x)-iB i {x). (3.26) 
Using (3.23) and (3.24) we also find 

2tr]T(^B l ) 2 = 4|(0„ -tgCM 2 + 4\(d, -igCM 2 + 4|d M 3 | 2 . (3.27) 

i 

Since 0i and 02 couple to C M in the same way and 3 does not couple to C M at all, we can 
choose 0i = 2 = = .B — iZ?, and 3 = i? 3 , so that 

2tr J2( V ^) 2 = 8j(^- ^ M )0| 2 + 4(^5 3 ) 2 . (3.28) 

i 

At large distances from the center of the flux tube, using cylindrical coordinates p, 9, z 
we have the boundary conditions: 

C -> - — , -> S e^, 5 3 -> S , as p -> oo. (3.29) 

The non-vanishing of i?o produces a color monopole current confining the electric color flux. 
The line integral of the dual potential around a large loop surrounding the z axis measures 
the electric color flux, just as the corresponding line integral of the ordinary vector potential 
measures the magnetic flux in a superconducting vortex. Since the dual potential is along 
a single direction in color space path ordering is unnecessary and the boundary condition 
(3.29) for C gives 

e -*9$ lmv e-<£ =e 2nir = e Mi) > (3 30) 

which manifests the unit of Z 3 flux in the tube. (A continuous deformation in SU (3) of our 
particular solution into a non-Abelian configuration will leave unchanged the path ordered 
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;a{h (i0i 



-Bl) 



+ (102 



-Bl) 



+ (103 



-Bl) 



integral Pexp(—ig J C-di) = e 2 ""^). The energy per unit length in this flux tube is the string 
tension a. The quantity g 2 / A plays the role of a Landau-Ginzburg parameter. Its value can 
be obtained by relating the difference between the energy density at a large distance from 
the flux tube and the energy density at its center to the gluon condensate. This procedure 
gives g 2 /X = 5 (which is near the border between type I and type II superconductors). We 
get from the numerical integration of the static field equations^ 17 ' 

(7^1.1(2452). (3.31) 

2 

We are left with two free parameters in C e ff, which we take to be a s = |^ = p and the 
string tension a. 

To couple C M to a qq pair separated by a finite distance we must represent quark sources 
by a Dirac string tensor G^ u . We choose the dual potential to have the same color structure 
(3.23) as the flux tube solution. Then G£, must also be proportional to the hypercharge 
matrix 

= YG^, (3.32) 

where G^L is given by eq. ()3.1()j) . so that one unit of Z 3 flux flows along the Dirac string 
connecting the quark and anti-quark. We then couple quarks by replacing G Mi , in C e ff ()3.19j) 
by 

= YG^ U , (3.33) 

where 

G^v = dfj,C v — dyC^ + G^ u . (3.34) 

Inserting (3.33) into ()3.19j) and using eq. ()3.28j) then yields the Lagrangian C e ff(G^ u ) 
coupling dual potentials to classical quark sources moving along trajectories z\(t) and z 2 {t): 

C CII (GJ - -3 4 + 2 + — W ' (3 35) 
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where W is given by (|3.25|) with cj>i = 2 = 0, 03 = -Ba- 
it is useful, as in (3.9), to decompose into its color electric components D and color 
magnetic components H. Similarly we decompose G^ u into its polarization components D$ 
and its magnetization components Hg'- 

D s = 2 e kmnG Smn , Hg = G$ k . (3.36) 

Then eq. (3.34) becomes 

D = -V x C + D s , H = -VC ~ — + H s . (3.37) 

The Lagrangian density £ e ff(G^) (3.35) can then be written as the sum of an "electric" 
part Co and a "magnetic" part C 2 . i-e., 

£eff(G„ v ) = C + £2, (3.38) 



where 



Co = - {^D 2 + 4|(V + igC)<P\ 2 + 2(V5 3 ) 2 + , (3.39) 



and 



C 2 = ^H 2 + A\(d - igC )<P\ 2 + 2(d B 3 ) 2 , (3.40) 

and all terms involving time derivatives appear only in £ 2 - 

We denote by W e ff(T) the Wilson loop of the dual theory, i.e., 

fVC ll V<j>VB 3 e'} dI " : ->>< G !~-' n+c <'rt 
The functional integral W e ff(T) eq. (3.41) determines in the dual theory the same physi- 
cal quantity as W(r) in Yang-Mills theory, namely the action for a quark- ant iquark pair 
moving along classical trajectories. The coupling in C e ff{Gf lv ) of dual potentials to Dirac 
strings plays the role in the expression ()3.41|) for W e ff(T) of the explicit Wilson loop integral 

e- !e i dlM " W in the expression (2.2) for W(V)^ 

4 We emphasize the distinction between W e ff(T) and the Wilson loop of the dual theory defined as an 

average of e 9 § ^ dl '. This dual Wilson loop would describe the interaction of a monopole antimonopole pair. 
For large loops the dual Wilson loop satisfies a perimeter law in accordances with 't Hooft's observation. 
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The assumption that the dual theory describes the long distance qq interaction in Yang- 
Mills theory then takes the form: 

W(T) = W eff (T), for large loops T. (3.42) 

Large loops means that the size R of the loop is large compared to the inverse mass of the 
Higgs particle (monopole field) <fi. Furthermore since the dual theory is weakly coupled at 
large distances we can evaluate W e ff(T) via a semi classical expansion to which the classical 
configuration of dual potentials and monopoles gives the leading contribution. This then 
allows us to picture heavy quarks (or constituent quarks) as sources of a long distance 
classical field of dual gluons determining the heavy quark potential. Thus, in a certain sense 
the dual gluon fields mediate the heavy quark interaction just as the electromagnetic 
field mediates the electron positron interaction. 

Using the duality hypothesis, we replace W(T) by W e ff(T) in eqs. (|2.4jl -(2.6) to obtain 
expressions for V (R) and Vyr> in the dual theory as the zero order and quadratic terms in 
the expansion of ilogW e ff(T) for small velocities 1\ and z 2 , i.e., the interaction Lagrangian 
Lj, calculated in the dual theory, is obtained from the equation 

i\ogW eff {T) = - f U dtLi{z u z 2 ,'z u 'z 2 ). (3.43) 

Remark 

There has been a recent revival of interest in the role of electric magnetic duality due 
to the work of SeibergJ 7 ' Seiberg and Witten' 8 ' and others on super symmetric non-Abelian 
gauge theories. Seiberg^ considered SU (N c ) gauge theory with Nf flavors of massless quarks. 
Although he did not exhibit an explicit duality transformation he inferred the complete 
structure of the magnetic gauge group and hence the associated massless particle content 
of the dual Lagrangian. For a certain range of Nf the dual theory is weakly coupled at 
large distances and hence the low energy spectrum of the theory consists just of the massless 
particles of the dual Lagrangian. Since this dual "magnetic" Lagrangian describes the same 
low energy physics as the original Lagrangian, the particle spectrum, mirroring the magnetic 
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gauge group, must appear as composites of the original "electric" gauge degrees of freedom. 
For Nf = N c + 1 the dual gauge group is completely broken, the associated dual gauge 
bosons become massive and the quarks of the original theory are confined. 

There are obvious differences between Seiberg's example, where the number of massless 
fermions plays an essential role, and the example of Yang-Mills theory where neither the 
original theory nor the proposed dual Lagrangian £ e ff contains fermions. Here confine- 
ment manifests itself via the development of a linear potential between heavy quark sources, 
whereas in the supersymmetric models confinement manifests itself via the realization of the 
hadron spectrum as composites of the original quark variables. In the supersymmetric model 
these hadrons are massless and as usual the production of these particles prevents the devel- 
opment of a linear potential. However, all the gauge bosons of the dual theory are massive 
and the coupling of the pure gauge sector to quark sources would produce a long distance 
linear potential between these sources. The common feature of Seiberg's supersymmetric 
model, where duality is "inferred", and Yang-Mills theory, where duality is assumed, is that 
in both cases the dual gluons receive mass via a Higgs mechanism which is the essential 
element of the dual superconductor mechanism. 



We now express the spin dependent heavy quark potential Vsd (2.8) in terms of quantities 
of the dual theory. As a first step we find relations of matrix elements of the dual field tensor 
Gfiv t° variations of W e ff(T) which are analogous to eq. (2.20) relating ((F^)) to variations 
in W(r). Consider the variation in W e ff(T) produced by the change 



4 The Potential V q q in the Dual Theory 





From eq. (|H.41jl we find that the corresponding variation 8W e ff(T) is given by 
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where 

«f(r ARK - 1 VC,V<pVB 3 e*f d ^ff^ +c ^) f{c ^ 0; Bs) 

<Pi B Z)l)eff = ~TTTr inS \>r ^ " K^' 6 ) 

JVC^V0VB 3 e l J M£*ff(G%„)+c G F) 
Using (|3.10|) to express the variation of G^ u in terms of the variation of the world sheet 
y M (cr, r), we obtain 

SG S 

dx-^(x)((G^(x))) ef , = 



(4.4) 



2 

The right hand-side of eq. ()4.4j) arises from varying the boundary of the Dirac sheet. The 
variation of the interior of the sheet produces a contribution proportional to the monopole 
current jff . 

3 MON {x) = d » Gfiv{x y (4 .5) 

(See eq. A. 52 of Reference 1 for details). This gives no additional contribution to eq. (|4.4|) 
since the monopole current must vanish on the Dirac sheet, so that no monopole can pass 
through the Dirac string connecting the charged particles. This latter assertion is just the 
dual of Dirac's condition for the consistency of a theory containing both electric charges and 
monopolest 14 ' . 

Defining dz x = dr^r, we can then write eq. f)4.4j) as 

5G S (x) 

dx^l{{cr{x))) ttSf = 



- e I (5z?dz^{(G Xa (z 1 ))) eff - 6z%dz*((G Xa (z 2 ))) eff ), (4.6) 



where 



G^(x) = l -e, vXa G x °(x). (4.7) 



Choosing a variation which vanishes on the curve T 2 , we obtain 



tfitogW e// (T) = --e J ^-^.((G^hff, (4.f 



where 8S tJ,v {zx) is given by eq. (j2.19|) . Eqs. (4.2) and (4.8) can be written as 

^W eff {T) = A e(( G (*)» ff = — — b x ^ l0g ^ (F) (4 9) 

which is the dual theory analogue of eq. ()2.20|) . The duality assumption ()3.42|) then gives a 
corresponding relation between matrix elements: 

((F^i))) = l((G^i)))ef f . (4.10) 

Eq. (4.10) gives a correspondence between local quantities in Yang-Mills theory and in 
the dual theory. The utility of electric-magnetic duality is that for large loops semi-classical 
configurations dominate the right hand side of eq. ()4.1()j) in contrast to the rapidly fluctuating 
configurations of Yang-Mills potential which contribute to the left hand side. Eq. ()4.10|) 
breaks up into its electric and magnetic components: 

- ((F mn )) = \e mnl {{H%„, (4.11) 

{{F M )) = \{{D t )) eff . (4.12) 

or equivalently, 

= \{{H e )) eff , (4.13) 

{{F m n)) = -e mn k{{D k )) eff . (4.14) 

Using eqs. (4.13) and (4.14) in eq. ()2.12|) gives the following expression for V^ S AG in the 
dual theory: 

V L M S AG = ~ £ ■ («#&)».// - *i x <0C*)».//), (4-15) 

where e\ — e and e 2 = — e. Note that {{H)) e jf — Vj x {{D)) e ff is the magnetic field at the 
jth quark in the comoving Lorentz frame, V^g AG the magnetic interaction of this field with 
a quark having a g factor 2. The fact that heavy quarks interact with a Dirac magnetic 
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moment is a consequence of the — expansion ^ for the qq Green's function upon which this 
analysis is based. 

To evaluate Vrhomas (2-13) we note from (j4.12|) that 



e 



2m? 
4 e 



3 2mf 



dtSffr x ((D( Zl ))) eff , (4.16) 



and obtain 



VThomas = ~ E \~ Sj ' (*& X 0(^)» 9// ). (4.17) 

The expression (4.17) is the contribution to the potential due to the precession of the axis of 
the comoving frame. In appendix A it is shown that ()4.17j) can be written in the usual form 

V Thomas = ■ * x R - -L^ a ■ v 2 x R. (4.18) 

Eq. ()4.18|) is essentially a kinematic relation and is independent of the dynamics of Yang- 
Mills theory. On the other hand V^| j4G , eq. (|4.15|) . depends upon the dynamics and cannot 
be expressed solely in terms of the central potential. 

To express V$s (2.17) in terms of quantities involving the dual theory we need the fol- 
lowing: 

zemiF^F^)}} - «^o(^)»«^o(4)»} 

_ A 6{(H k ( Zj ))) eff 

3 6 SHstty) ' 1 yj 



To obtain ()4.19|) we use eqs. (|2.21j) and (4.13) and the equation 

5((H k (Zj)))eff = eimnHHkjZjWeff 

2 5S™{ Zj ,) ■ 
(Compare eq. (4.9)). Using eq. (4.19) in (2.17), we obtain 



(4.20) 



/ MV SS = -i ± -f-T. ! it ! ( M . (4.21) 
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The factor multiplying SjSj, is symmetric in k and i and hence the terms in eq. (j4.21|) where 
j = f involve the combination 



q£ qk j_ qk q£ -i 

X k£ 



3 3 3 3 _ ~ X 

_ Q 



2 4 
Eq. (|4.21|) then becomes 



/ 



4^ e 2 f Mf ,J ((H k (z 3 ))) e ff 



3=1 3 "'3 "-"SKV/ 



3\mim 2 /Jri Jr 2 dHse{z' 2 ) 
The first term in eq. ()4.22j) is a spin independent velocity independent contribution to the 
potential proportional to inverse square of the quark masses while the second term in eq. 
(I4.22|) yields a spin-spin interaction of the expected structure. 
Finally, let us come to Voarwin (2-15) and note that 

((D"F u ,(z j ))) = d"{{F u ,(z j ))). (4.23) 

The derivative of the Wilson loop occuring in the definition (2.10) of ((F Ufl (x))) yields the 
Yang-Mills potential A v appearing in D"F vli . Using (4.10) we obtain! 

/ V Darmn dt = ^2 / dx^d u ((G^(x))) eff . (4.24) 

3 J •> 

For an alternative expression for Voarwin based on (eq. A. 7) of Ref. [4] see eq. B.3 of 
Appendix B. 

5 The Classical Approximation for Vo(R) and Vvd 

In the classical approximation eq. (3.41) becomes 



ilogW eff = - I dx£ eff (G s ), (5.1) 



5 Notice that ((F p „(z))} depends not only on the point z but on the entire Wilson loop. So in or- 
der for eq. (4.23) to make sense one has to use the appropriate definition of derivative. Given a 
functional 3>[ 7ab ] of the curve jab with ends a and b, under general regularity condition the varia- 
tion of $ consequent to an infinitesimal modification of the curve 7 — > 7 + Sj can be expressed as 
the sum of various terms proportional respectively to 6a, 5b and to the elements 5S pa (x) of the sur- 
face swept by the curve. Then the derivatives d / da p , d j db p and 5/5S f " 7 (x) are defined by the equation 
<5$ = d$/da p Sa p + d$/db p 5b p + J^SS f " 7 (x)S^/SS pa (x). In our case this would amount to put naively 

d/dzPPf(£dx»A ll (x)) = -Pf(J b z dx^A p (x))A p (z) and d/dz» f* dx^A^x) = A p (z)Pf'(f* dx^A^x)). 
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where C e ff(G^ u ) (|3.35|) is evaluated at the solution of the classical equations of motion: 

d a {d a Cp - dpC a ) = -d a G s af3 + jf ON , (5.2) 

(S M -^C7 M )V = ~^, (5-3) 

and 



where the monopole current jff ON is 



jf on = _ 3 ^[0*(^ _ igcjj - 0(«9 M + igC,)<f>*}. (5.5) 

As a result of the classical approximation all quantities in brackets are replaced by their 
classical values 

«GV(z)»e// = GW(a;). (5.6) 

The electric and magnetic components of eq. (5.6) are 

((D(x))) eff = D(x), ((H(x))) eff = H{x), (5.7) 

where D and H are the color electric and magnetic fields respectively given in terms of the 
dual potentials by eq. (3.37). 

We choose the Dirac string to be a straight line L connecting the quarks. As x approaches 
the string, (f>(x) — > 0, C^(x) — > C®(x), satisfying eq. (j3.17|) . As x — > oo, <p(x) — > B , C M (x) — > 
0, in contrast with the large distance boundary condition for the infinite flux tube. We can 
then choose <p(x) to be real so that 

cf>(x) = B(x), j™ ON (x) = -6g 2 C,B 2 . (5.8) 

Consider first the case of static quarks, v± = = 0. Then the scalar potential Co and 
the color magnetic field H vanish, and C e fj reduces to £ ec L- (3.39) which yields the static 
potential: 



V (R) = - I dxC , (5.9) 
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where £ is evaluated at the static solution of eqs. (5.2) - (5.4), which have the following 
form in this case: 

-V x (V x C) -Qg 2 B 2 C = -V x D s , (5.10) 

(-V 2 + g 2 C 2 )B = -y£(25£ 2 + 1B\ - 32E 2 ), (5.11) 

and 

4A 

- V 2 5 3 = -y5 3 (75 2 + - 16B 2 ), (5.12) 

where we have used the explicit form, eq. ()3.25|) . of W. 
To solve eq. (|5.1()jl it is convenient to write 

C = C D + c, (5.13) 

where C D is the Dirac potential satisfying the static form of eq. ()3.17j) namely: 

- V x (V x C D ) = -V x D s , (5.14) 

with Ds given by eq. (3.36). In cylindrical coordinates with the z axis along the line joining 
the two quarks at z = ±-R/2, eq. (3.10) gives 

D s = ee z {6(z - R/2) - 6{z + R/2)}5(x)5(y), (5.15) 

which describes the polarization vector for a line of dipoles. The solution of eq. ()5.14|) is 

C D = e^C D , (5.16) 

where 

c d_ e \ z-R/2 (z + R/2) ■ 

\ ^p 2 + (z- R/2) 2 .<Jp 2 + (z + R/2) 2 

Then 



c = e^c , (5.1J 
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and eq. ()5.1U|) becomes the following equation for c: 

(V 2 - 6g 2 B 2 )c = Qg 2 B 2 C D , (5.19) 

where 

^•^IKH+i- < 5 - 2o) 

Eqs. (|5.11|) . ()5.12|) and ()5.19|) are three nonlinear equations for the static configuration 
c, B, and B% with boundary conditions: c — > —C D ,B — > 5q, B3 — > i?o at large distances; 
c — ► 0, B — > for x on L. These equations have been solved in refJ 19 ' with the following 
results: 

The monopole current in eq. (|5.1()jl screens the color electric field produced by the quark 
sources so that as the quark anti-quark separation increases the lines of D are compressed 
from their Coulomb like behavior at small R to form a flux tube, and thus Vq(R) — > aR at 
large R. Both this small R and this large R behavior of the potential have their common 
origin in the evolving distribution of the flux of D whose divergence is fixed by the color 
electric charge of the quarks (V • D — V ■ D$) and whose curl is determined by the monopole 
current. Thus, the dual theory already in the classical approximation gives a potential which 
evolves smoothly from the large R confinement region to the short distance perturbative 
domain. This shows how the dual theory realizes the Mandelstam 't Hooft mechanism. 
It does not describe QCD at shorter distances where radiative corrections giving rise to 
asymptotic freedom and a running coupling constant are important. 

To calculate the terms in i log W e ff which are quadratic in the quark velocities we solve the 
field equations for moving quarks. To first order in the velocities the static field distributions 
follow the quark motion adiabatically. The time dependence of C, B and B3 then results 
from the explicit time dependence of R. Furthermore, since / dxCo generates the static field 
equations, it is stationary about the solution to these equations and remains unchanged to 
second order in the velocities. The velocity dependence in the potential then comes from 
the "magnetic" contribution C 2 which depends quadratically upon d C, d B 3 , and C , all of 

24 



which are first order in the velocities. The scalar potential C satisfies the equation, obtained 
from the time component of eq. (|5.2j) . 

V 2 C - 6g 2 B 2 C = V-H S , (5.21) 

valid to first order in the velocities. With the Higgs field B(x) already determined by the 
static equations, eq. ()5.21|) is a linear equation for the scalar potential, giving C to first 
order in the velocity. The velocity dependent potential Vvd is then given by 

V VD = -JdxC 2 , (5.22) 

representing the magnetic color energy due to the fields following the moving quarks. 

For small R the potential Vvd approaches the velocity dependent part of the Darwin 
potential (3.7) (multiplied by the color factor 4/3) because for small R the color magnetic field 
H(x) becomes the ordinary Biot-Savart magnetic field. As R increases the color magnetic 
field lines are compressed so that for large separation Vvd becomes linear in R. As an example 
consider the case in which two equal mass quarks move in a circular orbit of frequency uj. 
Then V\ = —v^ = i£2 ^ : , an d Vvd reduces to 

V VD = -\l{R)u\ (5.23) 

which defines the momentum of inertia I(R) of the rotating flux tube distribution. Eq. 
(J5.22j) evaluated for this configuration of moving quarks then determines I(R). For Large R 
we find' 1 ' 

lim I(R) = \{AR)R\ (5.24) 

it— >oo Z 

where 

A ~ .21cr, (5.25) 

determined numerically from eq. ([5.22)1 . By comparison we note that the moment of inertia 
V of an infinitely thin flux tube of length R is 

I'(R) = ^(A'R)R 2 , (5.26) 
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with 

A! = a/6. (5.27) 

The comparison of eq. (5.27) describing an infinitely thin flux tube with eq. (5.25) gives 
a quantitative estimate of the increase of the moment of inertia of the flux tube due to its 
finite thickness. 

We now compare these results for V + Vvd = —Li of the dual theory with the "modified 
area law" model' 9 !, eq. (2.22). In the dual theory i\ogW(T) is replaced by i log W e ff(T), 
given in the classical approximation by eq. (5.1). This gives in the limit of short distances 
the perturbative expression eq. (3.7) so that the short distance limit of the dual theory is 
the short range component i\ogW SR (T). The long distance limit of i\ogW e ff(T) is fixed by 
the values of a and A. Replacing A by A' in this limit yields % log W LR (T). This shows that 
i\ogW LR {T) describes a zero width flux tube. Aside from this difference we see that the 
"two components" of eq. (2.22) arise as two limits of a single classical solution describing 
the evolution of the potential produced by compression of the field lines with increasing R. 

As the simplest example of the implications of Vvd, we add relativistic kinetic energy 
terms to — (V + V V d) to obtain a classical Lagrangian, and calculate classically the energy 
and angular momentum of qq circular orbits, which are those which have the largest angular 
momentum J for a given energy. We find' 20 ' a Regge trajectory J as a function of E 2 which 
for large E 2 becomes linear with slope a' = J/E 2 = 1/8<t(1 — A/a). Then (5.25) gives 
a' 1/6.3(7, which is close to the string model relation a' = This comparison shows 
how at the classical level a string model emerges when the velocity dependence of the qq 
potential is included. The fact that the difference between the two expressions for a' is small 
indicates that the infinity narrow string may be a good approximation to the finite width 
flux tube forms between the qq pair. 

To summarize: 

1) The potential V (R) is determined by eqs.(5.9) and (3.39) evaluated at the static 
solution. 
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2) The potential Vvd is given by eqs. (5.22) and (3.40) evaluated at the solution of the 
classical equations to first order in the velocity squared. The resultant integrals have been 
calculated numerically W and determine four functions V + (R), V_(R), Vl(R) and V\\(R) which 
specify uniquely the terms in the potential proportional to the velocity squared. Explicit 
expressions for these functions are given in reference 1. 
Remarks 

1. In the absence of quark sources = 0), £ e // describes a system of massive dual gluons 
and monopoles. Because of the dual Higgs mechanism there are no unwanted massless 
particles in the spectrum. The massive particles of the dual theory cannot be identified with 
the massive particles of Yang-Mills theory, since the dual theory just describes the low energy 
spectrum. These masses determine rather the scale Rft — -h above which the dual theory 
should describe the qq interaction. Since a quark anti-quark pair moving in an orbit of radius 
R can only radiate a particle of mass M if ^ > M, in the domain R > jj where the dual 
theory describes Yang-Mills theory no dual gluons or monopoles are emitted. The glueballs 
of Yang-Mills theory, on the otherhand, are described by closed loops of color flux, obtained 
by coupling the dual potentials to closed Dirac strings and finding the corresponding static 
solution of the field equations of the dual theory. 

2. The Lagrangian density £ e // (3.35) describes the coupling of the Dirac string to Abelian 
configurations of dual potentials, and the functional integral (3.41) for W e ff(T) is restricted 
to such configurations. The external qq pair has in effect selected out a particular sector of 
the dual theory relevant to the qq potential. As a consequence the resulting potential should 
not be very sensitive to the details of the dual gauge group. 

3. The Dirac string in the classical solution was a straight line connecting the qq pair. This 



gave the configuration having the minimum field energy 



The flux tube corresponding to a 



given string position is concentrated in the neighborhood of that string since the monopole 



6 The Dirac string of the dual theory, in contrast to that of electrodynamics, is physical. The vanishing of 
the Higgs field on the string produces a vortex and an associated flux tube containing energy. This vortex 
can not be removed by a gauge transformation since such a transformation leaves the magnitude of the Higgs 
field unchanged. 
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current vanishes there. To evaluate the contributions to the potential arising from fluctua- 
tions of the shape and length of the flux tube' 21 !, we must integrate over field configurations 
generated by all Dirac strings connecting the qq pair. This amounts to doing a functional 
integral over all Dirac polarization tensors G~ a/ (x). Similar integrals have recently been car- 
ried out by Akhmedov, et al.J 10 ' in a somewhat different context. The functional integral 
over G^ u (x) is replaced by a functional integral over corresponding world sheets y(a, r), mul- 
tiplied by an appropriate Jacobian. As a result they obtain ' 10 1 an effective string theory free 
from the conformal anomaly' 22 '. Such techniques when applied in the context of the dual 
theory should lead to a corresponding effective string theory. 

6 The Classical Approximation for Vsd 

In this section we evaluate the expression for Vsd given in Section four using the classical 
solutions to the dual theory described in Section five. We consider separately the four 
contributions to Vsd (See eq. (12. 8j) ): 

(1) V Thomas : Eq. fnfl -(4.18) with V (R) determined by eq. ((Ojl . 

(2) Vlg AG : Eq. P~T5j) with ((D)) eff and ((H)) eff replaced by their classical values D and 
H, namely 

V L M S AG = ~ £ ■ 0(z 3 ) - v 3 x D(z,)\ (6.1) 

with H — vj x D calculated to first order in the velocity. To this order the static field 
configurations follow the motion of the quarks adiabatically and we find from eq. (3.37) 



H(zj) - Vj x D(zj) = -V(C (f) - C(x) ■ v(x)) 

where 



(6.2) 



and 



-V-A "I + v 2 _ _ 

v[x) = — hwxi, (6.3) 



Q = — ^L. (6.4) 
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In eqs. (|6.2H6.4|) . 2i±a j s the instantaneous velocity of the origin of the coordinates which we 
have chosen as the midpoint of the line L connecting the qq pair and u is the instantaneous 
angular velocity of L. (The motion of the qq pair along L does not contribute to eq. (|6.2Jl ). 

We can understand the result eq. ([6.2)1 . as follows: The left hand side is the color magnetic 
field at the position of the jth quark in the Lorentz system in which it is instantaneously 
at rest. The magnetic field in this co moving system is determined by the gradient of the 
corresponding dual scalar potential, namely Co — C ■ v. Indeed (6.2) remains valid beyond 
the classical approximation with the replacement Co — > ((Co)) e //,C — > ((C)) e //. 

Choosing R to lie along the z axis and using eqs. (5.21) and (5.10) for Cq and C we find: 

C Q -C-v = e 4> - C ^C-{z,p), (6.5) 
where p, <fr, z are cylindrical coordinates, and 

C-(z,p)=C D (z,p)+c^z,p), (6.6) 

where 

C D (z,p) = I ; 1 - 1 1 , (6.7) 

and where C-(z,p) satisfies the equation 

(V 2 - 6g 2 B 2 )c. = 6g 2 B 2 C D . (6.8) 

The solution of the linear integral equation ()6.8|) for c_ determines, via eqs. ()6.2|) and 
(6.5) the non-perturbative part of the color magnetic field in the comoving Lorentz system. 
From eqs. (6.6)-(6.8) it follows that for any fixed value of z and p this field vanishes like 
for large qq separation. The vanishing of this field at large R is in accordance with the 
observation of Buchmullert 23 ! that in a flux tube picture the color field in the comoving frame 
should be purely electric. However, for any finite value of the qq separation there is a color 
magnetic field in this system, and eqs. (6.1)- ([6.8|) give 



v mag = YMl I ( SrjRxp,) _ S 2 -(Rxp 2 ) \ + ( S 2 -(Rxp 1 ) _ S v (Rxp 2 ) 
LS R I I m 2 m 2 ) I m 1 m 2 mim 2 



29 



(6.9) 



where 



(6.10) 



p=0 

2 = fl/2 



The first term in (jfj.lOjl is the perturbative contribution to V 2 (R) arising from (7? and the 
second term is the non-perturbative part which behaves like for large R and which would 
not be present in the simple flux tube picture of Buchmuller. 



Finally adding Vl 4 s AG to Vrhomas gives the complete expression for the spin orbit coupling 



V LS , 
V ls 

where 



i dv 1 2 v;(R) 

R dR R 



S v Rxp! S 2 -Rxp 2 



2m\ 



2ml 



R 



S 2 -Rxp 1 S 1 -Rxp 2 



m 1 m 2 



m 1 m 2 



V[{R) = V'(R) - 



dVv 
dR' 



(6.11) 



(6.12) 



Eq. 1)6.11)1 expresses the spin orbit potential in terms of the central potential and a single 
independent function V 2 \R) determined by the dual scalar potential Cq—C-v in the comoving 
frame. This result for Vls satisfies identically the constraints of Lorentz invariance (6.12) 
(The Gromes Relations). ' 24 ' Furthermore, since V 2 \R) — *> for large R, we have 

dV 



lim V{{R) 

ft— >oo 



dR 



-CT, 



(6.13) 



which is the value given by the flux tube model for all R. 



(3) Vss' Eq. (4.22) with <5 ^^^ e ^ ff replaced by ■ Since, to first order in the velocity, 

C is determined by D s alone (see eq. ()5.10j0 the ^ term in H does not contribute to its 
variational derivative with respect to H$ and eq. (3.37) gives 



SH(zj) 



SH k {x) 
$H se (x> 



S k£ S(x - x)5{t - f) - V fe 



SC (x) 
5H se (x< 



(6.14) 



The quantity |S&- in turn satisfies the equation obtained by taking the variational derivative 
of eq. (|5.21j) with respect to H$, namely 



(V" - Qg 2 B 2 ^ 



SCq(x) 
'SHsAx' 



V e 6(x- x)5{t-t'). 



(6.15) 
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The double integral in eq. ()4.22j) then becomes a single integral over t of the static quantity 



. We emphasize that this simplification obtains only in the classical approximation 



&H(zj) 

we are now considering. 

Eqs. dSHD and (IfHKj) give 

8H k (2 



6 ke 8(x - x') + V k V £ G(x , x' ), (6.16) 



SHsi(x') 

where the Green's function G(x, 2') satisfies 

(-V 2 + 6g 2 B 2 (2))G{2, x') = 6(2 - x). (6.17) 

G(x, x') is the potential at 2 due to a point charge at 2 in presence of the monopole charge 
density j^ ON (5.8) carried by B(x). Since B(x) approaches its vacuum value B as 2 —* oo, G 
vanishes exponentially at large distances, i.e., 

-m B \x-x'\ 

G(2,2)^- - - (6.18) 
x 00 Aix\x — x'\ 

where 

o o „o 67T ~ 7T a , . 

m| = 6/S 2 = — 5 2 « - — , (6.19) 

and where we used the result, a « 2ABq, obtained from the energy per unit length of 
the static flux tube solution. Using a value a s = .37 obtained from fitting the cc and 66 
spectrum' 1 ' we obtain m# ~ 640 MeV. 

Separating off the Coulomb contribution to G we have 

G = -—r^ — =rr + G NP , (6.20) 

where G NP satisfies the equation 

(_ V 2 + 6g 2 B 2 )G NP = - % /f 2( ?, . (6.21) 
v 7 Att\x-x'\ v 7 

Inserting eqs. ()6.16|) and ()6.20j) into eq. (4.22) gives 

Vss = V°P + vV s m \ (6.22) 
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where 



VSt = \^—{ (Si ■ &)6& - z 2 ) + • V)(S 2 ■ V')G(x, Sf 



3 VTL\m 2 



X = Z\,X—Z2 ■ 



(6.23) 
(6.24) 



The potential Vg^ n is the same as previously obtained' 11 '. At small R it approaches the 
usual perturbative spin-spin interaction, and at long distances it is exponentially damped 

1 IttiP' 

due to screening by the monopole charge. The spin independent contribution V ss of V$s 
depends upon R via the dependence in eq. (|6.21j) of G NP on B. It was not included in 
ref.M. 

(4) V Darwin : Eq.(4.24) with {{G^)) eff replaced by G^ , namely 

4 ^ e 

To evaluate (|6.25J) we note from eqs. (3.11) and (|3.34J) that 



/ V Darwin dt = -\ J2 ^2 / dxfdTd^x). (6.25) 



d u G^(x)=j^x), (6.26) 

where j^x) is the quark anti-quark current. The monopole current does not contribute to 

d v G Vil and V Darwin becomes, 

4 y^/ dafjJx) = --Y 
4-8m?./r i M ; 3^ 8m, 

J J J 3 J 

Omitting self energy terms we insert p{z\) = —eS(zi — z 2 ), p(z 2 ) = e5(z 2 — z-\) into eq. (6.27) 
and obtain 



/ V Darwm dt = -- Y ^2 J r dfj^x) = E ^2 / dtp( Zj ). (6.27) 



Vbarwtn= 4" ( — + — ) ^1-^)- ( 6 -28) 

o \Tn 1 m 2 



In Appendix B we show that the alternate form (B.3) for VDarwin reduces in the classical 
approximation to the same expression (6.28). 

There are two then spin independent terms proportional to f 

^* + The first is 

(6.24). The second is VDarwin (6.28). 
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To summarize: In reference 4 the coefficient of + in the velocity dependent 
potential was written as: 



Vst 2 + VDarwin = J ( \ + \ ) V 2 (V (R) + K(#)) (6.29) 

8 V m{ mi 



which defines V a . Eqs. (6.24) and (6.28) give 



V 2 K = V 2 V NP (R) - |e 2 V ■ V'G NP {!Z,Z) 



(6.30) 



where V£ rp (R) is the non-perturbative part of the central potential so that V a is determined 
by the non-perturbative dynamics of Yang-Mills theory. The first term in (6.30) is the color 
electric contribution to V a and the second is the color magnetic contribution. 
The spin dependent potential is then given by: 

V 8D = V LS + Vg n + I (— 2 + -U V 2 (V (R) + V a (R)) (6.31) 

8 \mj 1712/ 

with V LS given by (6.10) and (6.11), by (6.23) and V a (R) by (6.30). 

It should be emphasized that these results along with those of section 5 do not account 
for quantum fluctuations about the classical solutions. To account for these fluctuations we 
must return to the more general eqs. (4.15), (4.22) and (4.24) and (4.3) which determine V q q 
in the dual theory independent of the classical approximation. 

7 Conclusion 

We have shown how the analysis of ref.t 3 ] of the heavy quark potential V q q in terms of Wilson 
loops W(T) leads to the expression for the long distance behavior of V qq in terms of an 
effective Wilson loop W e ff(T) calculated in a dual theory describing a dual superconductor. 
The coupling of the dual theory to heavy quarks is then uniquely specified with spin and 
relativistic effects accounted for unambiguously to order f A , the highest order for 

° J \ mass quark ) 1 ° 

which the concept of a potential makes sense. 

The calculation of W e ff(V) in the classical approximations leads to expressions for the 
various terms in V qq with clear physical interpretations. These results coincide for the most 
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part with the expressions for V q q given by a previous dual theory calculation^ in which the 
qq motion was treated semi-classically. The present treatment gives an additional contribu- 
tion, (6.29) to Vqq, and the contribution of "Thomas precession" now appears automatically 
whereas in the semi-classical treatment it has to be put in by hand. 

We can use the Wilson loop W e ff(T) to approximate W(T) for large loops, i.e., for 
R > jj, where M is either the mass of dual gluon C M or of the monopole field Bi (about 
500MeV). However, since the dual theory gives a heavy quark potential which approaches 
lowest order perturbation theory for small R, it should remain applicable down to distances 
where radiative corrections giving rise to a running coupling constant become important. 

Most significant is the fact that we have obtained an expression for the qq potential in 
the dual theory which makes no reference to the classical approximation . Furthermore since 
the formulae of reference ^ are obtained starting from a relativistic treatment of the qq QCD 
interaction, the results provide a direct connection of the dual theory to QCD which could 
lead to an understanding of the constituent quark model on a more fundamental level. 

As a final remark we note that the dual theory we propose is an SU(3) gauge theory, 
like the original Yang-Mills gauge theory. However, the coupling to quarks selected out only 
Abelian configurations of the dual potential. Therefore, our results for the qq interaction 
do not depend upon the details of the dual gauge group and should be regarded more 
as consequences of the general dual superconductor picture rather than of our particular 
realization of it. The essential feature of this picture is the description of long distance 
Yang-Mills theory by a dual gauge theory in which all particles become massive via a dual 
Higgs mechanism. 
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Appendix A 



Notice that 



Sj = \^r^P^j)))eff (Ai) 

can be interpreted as the acceleration of the jth quark so that eq. (4.17) can be rewritten 

1 2 

V T homas = -« S S 3 ' ^3 X ( A2 ) 

which is the usual expression obtained from semiclassical considerations. To express Vxhomas 
in terms of the derivative of the static potential we first note from eq. (4.9) that 

m ]f;f ] - >*)»„,. (A3, 

5Ds{x) 3 

Now, using the fact that 

V 1 D s (x) = -el5(x-z 1 ), (A4) 

— * o 

where Vi = we have 



o 



Uj y SD s (x)^ V ^S{X) ^ 

= —eR-J dx^((D(x))) e ff8(x — Zi) 
= -elR((D(z 1 )))effR- 

Now by symmetry, ((D(zi))) e ff evaluated at the position of a quark must lie along R. Hence, 

-±e({D{Z 1 ))) eff = R^. (A6) 
Eq. (Al) then gives a x = -^f|, so that eq. (A2) gives (4.18). 

Appendix B 

Here we begin with an alternate form for V Darwin where A a does not appear explicitly. ^ 
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dx»((D»F vlx {x)))= [ tf dtV 2 V 

J to 



-*e 2 / da? f dx' a [((F^)F^(x'))) - ((F„ u (x)))((F™(x>)))} = 



= [ f dtV*V -ie* f dt f dt'[{{F Gk {z 3 )F« k {z' 3 ))) - ((F ok ( Zj ))) ((F ok (z ;.)»]. 
Eq.(4.12) and the relation between 5S^ U and variations of G s give 



6 nvok,^ 4e 2 /■ r J{{D k ( Zj ))) eff 



(Bl) 



Then using eq.(2.21) with ^ 2 replaced by z\ and e by -e we obtain 

/ dtV Darwm = E / dfWo/SmJ - / * / ^A ^fn^n^ ' 

y J 3 ^ 8m j ^sfc(^) ( B3 ) 

which gives a second form for Voarwin- The classical approximation to (B3) is obtained by 
replacing 

EEMhn by **®-s(t-o. 

SDs^) SD s (zj) 

This yields the expression 



2 



^Darwin ^ ] 



1 _ 2 „, D v 4 e 2 



rW (i2) 



(B4) 



_8m] uv 7 3 8m^A?k(*i) 
Following the same reasoning that led to eq. (A5) we obtain 

? = V '"° W ' (B5) 

(There is no perturbative contribution to the left hand side of (B5).) The second term in 
(B4) then cancels the non-perturbative part of the first term. Eq. (B4) then becomes 



_L., 1 \t7*( 4 e 2 \ _ - z 2 ) ( 1 

which coincides with (6.28). 



V D „„, = _ + _ VM --— = e'^-^ ^ + -g , (B6) 
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t = t; 

Figure 1: Wilson loop for the quark anti-quark system. 



